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STRUCTURAL DESIGN METHOD AND RECORDING MEDIUM 



BACKGROUND OF THE INVENTION 
5 1. Field of the Invention 

The present invention relates to a method for designing a 
structure and particularly to a method for generating a model at 
an initial designing stage. 
2. Description of the Related Art 
10 Various methods for designing a structure have been proposed. 

«? For example, Japanese Patent Laid-Open Publication No. Hei 3-224063 
til discloses an optimal designing system in which an interface is 
01 provided between an analytical calculation section and modeling 
s section, changes of design parameters for a structure and 
f& remodeling are automatically performed based on the analytical 
2* result, and an optimal model which satisfies the design 
specification is output after design iteration. 

In this method, modification of the model makes the design 
specification satisfied after each iteration. However, this 
20 reference does not describe any technique for generating the basic 
original shape or configuration. Current structural design depends 
on a designer's experience or on his design notes, and no tool exists 
for generating the original shape (analysis model). Certainly, 
it is possible to re-design a configuration using the above method 
25 or a design support program, such as the AISI/CARS program for 
automobile structure member design once the original shape is 
generated. 



Finite element analysis (FEA) based on detailed mesh models 
is widely used in the analysis phase. However, much calculation 
time and the number of design changes are required to regenerate 
a model after each design change in the case where the original 
5 model is unreasonable. This is because the convergence and 
calculation time is determined from the completeness of the 
original model. 

SUMMARY OF THE INVENTION 
The present invention was conceived to solve the above- 

IQ mentioned problem and one of the objects of the present invention 
is to provide a method for designing a structure which can easily 

m and quickly generate a reasonable shape at an initial design stage. 

fjl In order to achieve the object mentioned above, one aspect 

s of the present invention provides a method for designing a structure 

W comprising the steps of dividing the design domain of the structure 
21 into a plurality of nodes, coupling adjacent nodes and coupling 
non-adjacent nodes by beam elements to set a design model, and 
optimizing the design model. 

According to another aspect of the present invention, a 

20 boundary condition is set to the design model and stiffness of beam 
elements is reconfigured in the optimization process. Here, it is 
preferable that the boundary condition includes an external force 
applied to the structure, a fixed condition, and a volumetric 
condition, and the optimization performs substantial removal of 

25 the beam elements. For a design taking panels into account, it 
is preferable to add the stiffness contribution of panels to the 
design domain prior to the optimization. 



In the present invention, the design domain of the structure 
is discretized by a plurality of nodes with the nodes connected 
with each other by respective beam elements. A modeling scheme 
called as a ground structure approach is well known. In this scheme, 
5 a design domain is discretized by a plurality of nodes with the 
adjacent nodes connected with each other by respective trusses. 
This approach is referred to, for example, " Optimization of 
Structural Topology Shape and Material", Martin Philip Bendsoe, 
Springer, 1995, However, in the present invention, in addition 
|LQ to the adjacent nodes being connected as in ground structure 
0! approach, non-adjacent nodes are also connected and, moreover, 
01 these connections are performed by beam elements, rather than truss 
y» elements. Connecting the non-adjacent nodes can increase the 
f number of variations in the generated shape. Furthermore, by using 
M a beam element for connection in the present invention, rather than 
r l truss element as in the related models, more detailed and exact 
optimal shapes satisfying the predetermined limitation condition 
can be generated. All the adjacent nodes must be connected, but 
not all the non-adjacent nodes must be connected by beam elements. 
20 However, the accuracy of the generated model increases as the number 
of connected non-adjacent nodes increases. 

As a method for optimizing a design model constructed using 
nodes and beam elements, for example, a optimal shape with maximum 
stiffness can be generated by applying a predetermined boundary 
25 condition to the design model and modifying the stiffness of beam 
elements. As a boundary condition, an external force and fixed 
condition can be applied, and optimization can be performed so that 
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the displacement under a constant volume constraint is minimized* 
The optimization is performed by adjusting the stiffness of beam 
elements, and a predetermined volume constraint can be satisfied 
by substantially removing the beam elements which do not 
substantially contribute to maximization of stiffness. 
Substantial removal includes cases where the beam elements are 
directly deleted from the data and also cases where the practical 
existence of the beam elements is negated by setting zero to a 
parameter of the beam element which affects existence, such as, 
for example, the cross sectional area. When data itself is deleted, 
the number of data which plays a role in the calculation changes, 
resulting in complexities in the calculation algorithm, but by 
substantially removing data, the complexity in the calculation 
algorithm can be inhibited. 

When designing a structure such as an automobile body or a 
house, panel elements contribute to stiffness in addition to the 
frame structure. Thus, by adding the stiffness contributions 
from panel elements to the nodes in the design model constructed 
from nodes and beam elements, optimization can be performed using 
the same algorithm in the case where no panel element is employed. 

It is preferable that the cross section of the beam element 
has a shape that is represented by at least 2 design variables. 

By setting the number of design variables to 2, elliptic and 
rectangular cross sections can be represented, and by setting the 
number to 3 or more, various complex cross sectional shapes can 
be represented, with such a configuration, at the optimization 
step, design of beam cross sections is also possible in addition 



to the ability to lay out the beams. 

In particular, it is preferable to determine dimensions of 
the cross section and angle of the main axis of the cross section 
with respect to the local coordinate system of a beam element during 
optimization process. 

Therefore, it is possible to determine the main axis 
direction (angle) of the cross section in addition to the dimension 
of the cross section. This results in increasing the number of 
degrees of freedom for design and in improving the performance of 
the product. In practice, most of the beam structures used for 
machine parts have a cross sectional main axis. Thus, the 
capability to determine the direction of the cross sectional main 
axis is very useful in machine designing. 

Another aspect of the present invention provides a computer 
readable medium for storing a program for designing a model for 
a structure. The program is characterized in that, when executed 
by a computer, it causes the computer to discretize the design domain 
of the structure by a plurality of nodes, generate a design model 
which is constructed by connecting the adjacent nodes and 
connecting the non-adjacent nodes by beam elements, and optimize 
the design model. 

It is also preferable that the program, when executed by a 
computer, causes the computer to perform the optimization by 
modifying the stiffness of the beam elements based on an input 
boundary condition, where the boundary condition includes an 
external force applied to the structure, a fixed condition, and 
a volume constraint. It is also preferable that the optimization 



includes substantial removal of the beam elements* 



BRIEF DESCRIPTION OF THE DRAWING 

Fig. 1 is a structural diagram of a design model according 
to a preferred embodiment of the present invention. 

Fig. 2 is a processing flowchart of the preferred embodiment 
of the present invention. 

Fig. 3 is a descriptive diagram showing the result of 
optimization of the design model depicted in Fig. 1. 

Figs. 4A, 4B, and 4C are diagrams showing specific examples 
of cross sectional shapes . 

Fig. 5 is a diagram showing a design model for an automobile 

body. 

Figs . 6A and 6B are diagrams explaining loads to be added 
to the automobile model depicted in Fig. 5. 

Fig. 7 is a diagram explaining a model for optimizing the 
automobile body depicted in Fig. 5. 

Fig. 8 is a descriptive diagram showing a connection section 
which is the object to be designed. 

Fig. 9 is a diagram showing a design model of the connection 
section depicted in Fig. 8. 

Fig. 10 is a diagram showing a model for optimizing the design 
model depicted in Fig. 9. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 
A preferred embodiment of the present invention will now be 
explained referring to the attached diagrams. In practice, the 



design method of the present invention can be performed on a computer 
by installing the processing program to a computer such as a personal 
computer and a workstation. Because the design model of the present 
invention uses elements based on structural mechanics , such as a 
5 beam element and a panel element, a high capacity memory is not 
necessary, and it is possible to, for example, use common 
spreadsheet software as a front end interface on a personal computer. 
The processing program can be recorded in any preferred medium which 
can hold information electromagnetically , optically, or chemically, 
10 such as a CD-ROM, DVD, a common hard drive, or a semiconductor memory. 
5; The installation of the processing program can be executed by, for 
pi example, storing the processing program on a CD-ROM and supplying 
IP the processing program from the CD-ROM to the hard drive of a 
£ computer. It is also possible to store the processing program in 
fi§ a hard drive or a ROM of a computer and use the computer as a computer 
j;f dedicated to designing structure models. 

Fig. 1 shows schematically a design model used to illustrate 
the preferred embodiment of the present invention. The design 
domain of an object to be designed or an object to be analyzed 
20 (hereinafter, the term "design" is assumed to include analysis) 
is descretized by a plurality of nodes 10, and adjacent nodes are 
connected by a beam element 12. In addition, non-adjacent nodes 
are also connected by beam elements 12 . In Fig. 1, all the adjacent 
nodes are connected by beam elements 12 and all the non-adjacent 
25 nodes are also connected by beam elements 12. Connecting all 
non-adjacent nodes by beam elements 12 is not necessary, but it 
is advantageous to have the number of beam elements 12 as large 



as possible in order to generate an optimal initial design shape. 
In the example of the present embodiment, since an optimal original 
shape is generated even though no shape (configuration) is provided 
at the initial stage, the number of varieties of shapes which can 
be generated increases as the number of beam elements increases; 
thus, it is preferable to connect all the non-adjacent nodes by 
beam elements 12. 

In this manner, a design model is set, and the design model 
is then optimized by setting boundary conditions* 

Fig. 2 is a processing flowchart for maximizing the stiffness 
of a structure using the design model depicted in Fig. 1 as a model. 
First, an initial value is input (step S101). Here, "inputting 
the initial value" refers to discretizing the design domain by a 
plurality of nodes 10, connecting the nodes 10 by beam elements 
12, and inputting of the boundary conditions. The boundary 
condition includes an external force load and a volume constraint. 
In the present example, we describe a case where the structure on 
the left side in Fig. 1 is constrained and a load is applied to 
the central portion 14 at the right side. 

After the initial value is input, stiffness matrix K of the 
entire structure is constructed (step S102) . In general, when the 
design model comprises a plurality of nodes 10 and beam elements 
12, the relationship between force and displacement of nodes can 
be represented by f == Ku, where f is the force vector, u is the 
displacement vector, and K is the stiffness matrix. K in the 
equation for defining the relationship between the force and 
displacement is the stiffness matrix of the structure, and can be 
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represented by superposition of elemental stiffness matrices K e 
represented by the following equation, 



(Eq 1) 
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where A is the cross sectional area of a beam element, E is the 
Young's modulus, G is the shear modulus, 1 is the length of a beam 
element, i^, l az , and J are second moments of inertia. 
Specifically, Iyy, I 3Z , and J are respectively second moments of 
inertia with respect to bending in y axis direction and z axis 
direction, and torsion of the beam around the x axis. In thhis 
matrix, "Symmetric" refers to a symmetrical matrix. From the above 
equation, it can be clearly understood that a beam element is 
different from a truss element. More limitations arise when 
modeling using a truss because a truss has stiffness with respect 
to only one axis whereas a beam element has stiffness with respect 
to six axes, in the present invention, by using a beam element, 
instead of a truss element, an optimal configuration which has a 
greater degree of freedom, but which still satisfies a 
predetermined limiting condition, can be generated. 

After stiffness matrix K is generated, the sensitivity of 
K with respect to a design variable is obtained (step S103). The 
simplest example would be a case where the cross section is assumed 
to be a solid circle and the design variable is set as the cross 
sectional area. In such a case, the direction of the main axis 
of the cross section is not determined. A design variable in this 
example is the normalized cross sectional area p (0 S P S 1) and 
a relationship, 

(Eq 2) 
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holds with respect to the cross sectional area A, where is the 
maximum value, and p ultimately becomes either 0 or 1. Normalized 
cross sectional area p is used in order to simplify the calculation. 
The differential value of K with respect to p must be determined 
5 in order to determine the sensitivity of the stiffness matrix K 
with respect to the design variable p. If the cross section of 
a beam element is assumed to have a diameter of d, the following 
relationships are obtained. 

1ft (Eq 3) 

HI j «/ =± d *- p2A ™ 
™ ^ 64 4jv 

Jj! < E( 3 4 ) 



32 2x 



liy Thus, relationships yield, 



(Eq 5) 
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(Eq 7) 

Bp Jt 

By using these equations, the sensitivity of the stiffness matrix 
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K with respect to the design variable p can be derived. 

After the sensitivity of the stiffness matrix K is obtained , 
an equation of equilibrium, f = Ku, is solved for u (step S104). 
Specifically, by multiplying the external force f by an inverse 
5 matrix K" 1 of the stiffness matrix K, u can be determined. After 
determining u f the mean compliance, a volume constraint, and 
respective sensitivities are calculated (step S105). The mean 
compliance is an objective function which is the objective for 
optimization. In the case of maximizing the stiffness, it is known 
If)! that stiffness can be maximized by designating (f -u) as an objective 
% function, and minimizing the objective function (f-u). Here, 
d% vector f is an external force given as a boundary condition and 
u is derived at step S104. Therefore, by using these values, (f- 
u) can be determined. The volume constraint can be set as ratio 
l^s of initial volume of the initial configuration. For example, if 
?J the initial volume of the design domain is set as 100%, a volume 
limit of 20% (i.e. reduce 100% to 20%) can be given. The volume 
limit can also be given as a weight (for example, limit 100 kg to 
20 kg) . The sensitivity of the mean compliance is obtained as the 
20 following relationship. 

(Eq 8) 

f = fu=u T Ku 

25 Therefore, the sensitivity of (f*u), that is, the differential value 
of (f # u) with respect to the design variable p, is, 
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(Eq 9) 

3/ T dK 
dp dp 

where T represents transpose of a matrix. Because the sensitivity 
of the stiffness matrix K is already determined in the process at 
step S103, by using the sensitivity, the sensitivity of (f*u) can 

be determined. The sensitivity of the volume constraint can also 
be determined by determining the differential of the volume with 
respect to the design variable p • 

After determining the mean compliance and volume constraint, 
and respective sensitivities, the design variable p is updated and 
optimized by an optimization method using these values (step S106) . 
As an optimization method, any of the known methods can be used 
such as SLP (sequential linear planning) or CONLIN (a type of 
sequential convex function approximation) . The optimization 
algorithm describes a scheme to solve the problem of minimizing 
the objective function, (f*u), and is an algorithm to determine 
p which minimizes (f*u) considering respective sensitivities of 
the mean compliance (f*u) and volume constraint with respect to 
p. Specifically, the cross sectional area of a beam element 12 
which has a high sensitivity is increased and the cross sectional 
area of a beam element 12 with a low sensitivity is reduced or the 
beam element itself is substantially removed by setting the cross 
sectional area to zero. By substantially removing the beam 

elements 12 with low sensitivity, the volume can be reduced and 
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the volume constraint can be satisfied. After the optimization , 
it is then judged whether or not the optimization converged {i.e. 
a minimum solution is obtained) (step S107), and if not, processes 
at steps S102 through S106 are repeated- In other words, the cross 
sectional area of the beam elements 12 are changed, value of (f- 
u) is examined, and if a minimum solution cannot be obtained, the 
cross sectional area of the beam elements 12 are again changed for 
optimization . 

Fig. 3 shows a shape optimized through the processes 
described above for the structure depicted in Fig. 1. Here, only 
the beam elements 12 with high sensitivity remain, now with their 
stiffness optimized. In other words, these beam elements remain 
with increased cross sectional area, whereas the beam elements 12 
with low sensitivity are substantially removed to yield a shape 
with maximized stiffness according to the predetermined boundary 
conditions (constraint condition, external force, and volume 
constraint). The thickness of the lines in the figure stands for 
the size of the cross sectional area of the beam elements 12, where 
thicker beam elements 12 have larger cross sectional areas. Such 
a shape can be used as an initial design shape, and more detailed 
analysis can then be performed by FEA when preferable. Even when 
FE analysis is performed, because the initial shape depicted in 
Fig. 3 is a mechanically reasonable shape, the analysis will 
converge quickly and the calculation time will be short. 

While in the above description a beam element 12 whose cross 
section is a solid circle with a cross sectional area of 7Td 2 /4 
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was chosen as a simple example, the invention is not limited to 
simple beam elements and can be preferably employed with more 
complex shaped beam elements 12. 

Consider a case where k design variables are set for one beam 
cross section, each design variable represented by x L (x = 1, 2, 
k) with respective maximum value of x^. Using design variables 
Pi, which are normalized values for design variables x ir the 
relationship is represented as follows. 

(Eq 10) 

P L satisfies the condition, 

(Eq 11) 
Osp.xl 

During optimization, the design variables are replaced by 
these p L to determine the optimal solution. Moreover, when the 
cross sectional shape is not symmetrical with respect to two main 
axis directions, a main axis direction must be determined to 
minimize the objective function. Here, the angle of the main axis 
direction with respect to the local coordinate system of a beam 
is represented by 6 to determine an optimal rotational angle. 

The mean compliance shown in the following equation is used 
similar to the above case as an objective function f . In other 
words, with the load vector f applied to the structure and 
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displacement vector u, the objective function can be represented 
as, 

(Eq 12) 
5 f = fu=u r Ku 

similar to the above equation 8. 

The objective function f can also be represented by, 

Ifk (Eq 13) 

0! This is because the objective function is a function of normalized 

f l values (design variables) /0 ± (i = 1, 2,3,,.., k) and a rotational 

iu angle 8 . 

Ill The optimization should be performed to satisfy the following 

O equation. 

(Eq 14) 
minimize / « / *u - u T Ku 

20 

The constraints are, 

(Eq 15) 

25 0<: pi si fori=l,-,* 

Ku-f 
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where V is the total volume of the structure and V u is a preset maximum 
volume constraint. Then, we set A, A lf and A 2 as Lagrange 
multipliers to define the following Lagrangian L. 



5 (Eq 16) 

£-/+A(^-^)+^(- A )+^( A -l) 



The KKT-conditions are, 



lg s (Eq 17) 

S dL _ % 

f dp, dp, [dp,) 



dL df ( dV \ , T dK . ( dV \ . , n 
— -— + A| \ +Aj =-« r « + A \-\ +Aj = 0 

dp, \dPi) 

dL df T BK 
-« — u = 0 



p dd d8 dd 

H a(v-v u )=o 

% A,(a-1)-0 

H AzO 
I! Aj a 0 

a: 0 

A -1*0 

20 

The optimal solution can be obtained by determining design 
variables p L which satisfy these equations. Here, any of the known 

methods can be used for design variables p L , similar to the above 
case for a cross sectional area as one design variable. For the 
25 design variable 0 , on the other hand, an optimal solution can be 
obtained by determining a solution that satisfies the equation 11. 
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in the optimization method mentioned above , it is necessary to 
determine sensitivities of the objective function f and of the 
limiting conditions with respect to the design variables p± m The 
sensitivity of the objective function can be determined from, 
similar to the above equation 9, 

(Eg 18) 

df T dK 
— u 

dp. B Pi 

where K is the stiffness matrix of the whole structure. This 
equation is similar to equation 9 except that a plurality of 
p s must be taken into account* That is, in this example , since 
the number of design variables is i f sensitivities for each design 
variables must be determined. 

Figs. 4A through 4D show four cross sectional shapes to which 
the present invention can be applied. Fig. 4A shows an elliptic 
cross section, Fig. 4B shows a rectangular cross section, Fig. 4C 
shows a hollow rectangular cross section, and Fig. 4D shows a C-type 
cross section. 

For example, in the case of an elliptic cross section depicted 
in Fig. 4A, the cross sectional area A of the beam, second moment 
of inertia I y in the y axis direction, second moment of inertia l z 
in the z axis direction, and polar moment of inertia J can be 
determined as 
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(Eg 19) 
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Using these cross sectional area A of the beam, second moment 
of inertia I y in the y axis direction, second moment of inertia I z 
in the z axis direction, and polar moment of inertia J , the stiffness 
matrix K Q of a beam element can be derived from R and K ep as follows . 
(Eq 20) 

Ke = R T KepR 



where K ep and R are defined by 



(Eq 21) 
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When two design variables are used such as the cases depicted 
l#l in Figs. 4A through 4D, the design variables x 2 and x 2 can be 
y! represented by, 



(Eq 22) 



using normalized design variables p x and p 2f similar to the earlier 
case. By determining respective sensitivities of A, I y/ I z and J 
with respect to p x and p 2 , respectively, sensitivities of the 

20 stiffness matrix of the beam element with respect to the design 
variables, dK/dp 1 and dK/ d p 2 , can be obtained, and thus, the 
sensitivity of the objective function f can be computed using 
equation 18. From the sensitivity, optimization calculation can 
be performed as above, to obtain an optimal solution. 

25 The update value of the rotational angle 6 of the beam element 

(angle of the cross sectional main axis with respect to the beam 
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element coordinate system) can be determined from an equation 
derived from the KKT-conditions . 



(Eg 23) 
r dK 



u =0 

dd e 



The partial differential equation of the objective function f with 
respect to the design variable 6 yields polynomial solutions of 
sinusoidal functions, and thus, the value of 0 can be determined 
from any known value calculation method. 

In other words, K & becomes, from the above equation 20, 
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(Eq 24) 
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where I lf I 2 , and I 12 are 
(Eq 25) 

It follows that the sensitivity of K a is 
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(Eq 26) 
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where I = and I„ are 



(Eq 27) 



[l y -l z )sin20 



39 60 



= f 7 „ _/ , )cos20 =^12. 
\y 2 3d 



If the displacement vector U a of the element is represented 

by 

m (Eq 28) 



Ue 



-{»U "ly «U *U ^ "2, " 2 y «2x ^ ^ ^ ^ 



equation 23 becomes 



(Eq 29) 
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^ -I z j{c t cos20 + C 2 sin2#'} = 0 
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where C 1 and C 2 are 
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(Eq 30) 

C l =-6u ly u u + 6u ly u 2z + 6u u u 2y -6u 2y u 2z 
+3Lu ly d ly + 3Lu ly 6 2y ~3Lu u 6 u -3Lu lz 6 2z 
-3Lu Zy 9 ly -3Lu 2y 9 2y + 3Lu 2z d iz + 3lu 2z 9 lz 
+2L%0 U + L%0 2y +L\d 2z +2L\ y d Z! 

2 2 2 2 

C 2 = 3u ly -6u ly u 2y -3u u + 3u 2y + 6u u u u -3u 2z 
+3Lu iz d iy -3Lu 2z 6 ly -3Lu 2y d lz + 3Lu lz 6 2f 
-3Lu 2z d 2y -3Lu 2y 6 2z +3Lu ly d i2 + 3Lu ly 6 J 
+L% 2 -L%'S 2y - -L\' 2 -L\ y \L% 6 2z +L% Z 2 

Therefore, the optimal rotational angle 0 can be determined 
by the following procedure. First, if I y equals I z , the cross 
section would be symmetrical with respect to two axes and no 
rotational angle can be determined. Thus, we assume the l y and l z 
are different. 

If C 2 is equal to zero, cos (2 0) will be equal to zero, and 
thus, 0 is equal to either -7T/4 or 7t/4. 

If, on the other hand, C 2 is not equal to zero, a solution 
can be obtained from an equation, tan (20) = -CjVC^. 

In either case, two or more solutions are determined. The 
sufficient condition for the solution is 



(Eq 31) 



d z L = JTl 8D 
dd 2 dd 



T^'^T^y ~ l z j{-Qsin20 +C 2 cos20 }>0 



By evaluating this equation, a solution can be determined. 

In this manner, a rotational angle 6 can be determined for 



an elliptical case. Thus, the main axis direction of the beam 
element can be determined and cross sectional design of the beam 
element can be performed in addition to the layout of the beam 
element in the optimization step. 

Fig. 5 shows an example of applying the design method of the 
embodiment on an automobile body. Fig. 5 is a design model in which 
the design domain of an automobile body is constructed from a 
plurality of nodes 10 and beam elements 12. 

Figs . 6A and 6B show an example where the design model depicted 
in Fig. 5 is given boundary conditions, of a constraint in the right 
side as a fixing condition and of a load at the left side as an 
external force. Fig. 6A shows a case where a load is applied from 
top left of the top view and Fig. 6B shows a case where a load is 
applied from bottom left of the top view. An optimal stiffness 
is obtained in which the displacement is minimized when these loads 
are separately added (i.e. in so-called multi-loading). 

Fig. 7 shows a result of the optimization on a design model 
depicted in Fig. 5 with a boundary condition depicted in Figs. 6A 
and 6B, following the processing flowchart depicted in Fig. 2. 
Similar to Fig. 3, a reasonable initial shape is obtained with beam 
elements 12 having low sensitivity to the stiffness substantially 
removed and only the beam elements 12 having high sensitivity 
remaining. 

The application of the design method according to the present 
invention is not limited to only an overall structure such as an 
automobile body or a skeleton structure of a building, and can also 
be applied to a partial structure such as a connection section of 



a body skeleton. 

Fig. 8 shows an example partial structure to which the design 
method of the present invention can be applied. An initial shape 
for a connection section 2 0 in an automobile body can be generated. 

Fig. 9 shows a design model of a connection section 20, and 
the design model is constructed from a plurality of nodes and beam 
elements connecting the nodes, similar to Figs. 1 and 5. 

Fig. 10 shows a result of an optimization of the design model 
depicted in Fig. 9 following the processing flowchart depicted in 
Fig. 2 under a predetermined boundary conditions . Unnecessary beam 
elements having low sensitivity are substantially removed and a 
shape with the stiffness (that is, the cross sectional area) of 
each beam element optimized is obtained. It is possible to further 
analyze in more detail by FEA using the shape depicted in Fig. 10 
as an initial shape. 

Because in an automobile body, for example, panels exist 
outside the design domain which add stiffness to the form within 
the design domain, it is preferable when designing an automobile 
body to optimize the form considering the stiffness contribution 
by the panels. The stiffness contributions by the panels can be 
added as a property of the nodes within the design domain, and by 
performing the optimization using the nodes with the stiffness 
added and beam elements, a form with optimized stiffness can be 
obtained, which takes the panels into consideration. 

In the example used to illustrate the preferred embodiment, 
stiffness maximization was described as an example of optimization. 
However, the present invention can similarly be applied to 
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optimization of eigenf requency (a form with maximized 
eigenfrequency ) or of flexibility. The optimization can be 
performed by determining the sensitivities of the beam elements 
with respect to the design variables and changing the design 
variables considering the sensitivities. 
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What is claimed is: 

1. A method for designing a structure comprising the steps 

of: 

discretizing the design domain of said structure by a 
plurality of nodes; 

creating a design model by connecting adjacent nodes by beam 
elements and by connecting non-adjacent nodes by beam elements; 
and 

optimizing said design model. 

2. A method according to claim 1, wherein, 

a boundary condition is given to said design model, and 
said optimization is performed by correcting the stiffness 
of said beam elements. 

3. A method according to claim 2, wherein, 

said boundary condition includes an external force applied 
to said structure, a fixed condition, and a volume constraint; and 

said optimization includes substantial removal of said beam 
elements . 

4. A method according to claim 1, wherein, 

the cross section of said beam element has a shape which can 
be represented by two or more design variables. 

5. A method according to claim 4, wherein, 

a dimension of said cross section and an angle of main axis 
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of said cross section with respect to a coordinate system of said 
beam elements are determined at said optimization step. 

6. A computer readable recording medium for recording a 
5 program for designing a model of a structure, characterized in that 
said program, when executed by a computer, causes said computer 
to, 

discretize the design domain of said structure by a plurality 
of nodes; 

|C) create a design model by connecting adjacent nodes by beam 

m elements and by connecting non-adjacent nodes by beam elements; 

Cm and 

optimize said design model. 

15 7. A recording medium according to claim 6, wherein, said 

ri program, when executed by a computer, causes said computer to 

perform said optimization by correcting the stiffness of said beam 

elements based on an input boundary condition. 

20 8. A recording medium according to claim 7, wherein, 

said boundary condition includes an external force applied 

to said structure, a fixed condition, and a volume constraint, and 
said optimization includes substantial removal of said beam 

elements . 

25 

9. A recording medium according to claim 6, wherein, 
the cross section of said beam element has a shape which can 
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be represented by two or more design variables. 

10. A recording medium according to claim 9, wherein 
a dimension of said cross section and angle of main axis of 
said cross section with respect to a coordinate system of said beam 
element is determined at said optimization step. 
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ABSTRACT OF THE DISCLOSURE 

The design domain of an object to be designed is discretized 
by a plurality of nodes (10), and adjacent nodes and non-adjacent 
5 nodes are connected by beam elements (12). Optimization is 
performed on a design model constructed from the nodes (10) and 
beam elements {12), to obtain an initial design shape. When 
maximizing stiffness , the sensitivity of the stiffness matrix is 
determined, the cross sectional area of the beam elements with high 
10 sensitivity is increased, and the cross sectional area of the beam 
2; elements with low sensitivity is either reduced or set to zero so 
m that low sensitivity beam elements are substantially removed. By 
fp constructing the design model using a plurality of beam elements, 
h a reasonable shape satisfying the boundary condition can be 
IP obtained. 
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